The electronic structure of the heavy fermion metal LiV 2 04 
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The electronic structure of the first reported heavy fermion compound without f-electrons L1V2O4 
was studied by an ab-initio calculation method. In the result of the trigonal splitting and d-d 
Coulomb interaction one electron of the d configuration of V ion is localized and the rest partially 
fills a relatively broad conduction band. The effective Anderson impurity model was solved by Non- 
Crossing-Approximation method, leading to an estimation for the single-site Kondo energy scale 
Tk- Then, we show how the so-called exhaustion phenomenon of Nozieres for the Kondo lattice 
leads to a remarkable decrease of the heavy-fermion (or coherence) energy scale T co h = Tk 2 /D (D 
is the typical bandwidth), comparable to the experimental result. 



Heavy fermion systems are characterized by a large ef- 
fective quasiparticle mass inferred from the strongly en- 
hanced electronic specific heat coefficient and spin sus- 
ceptibility at low temperatures Until recently this 
effect was observed only for f-electron compounds con- 
taining lanthanide or actinide atoms. The transition 
metal oxide compound LiV2 04 is the first reported heavy 
fermion system without f-electrons Q . 

L1V2O4 has a face-centered-cubic normal-spinel struc- 
ture. The formal valence of the V-ions is ]/ 3 - 5 + lead- 
ing to 1.5 electrons/ V in the 3d-band. The electronic 
specific heat coefficient 7(T) = C e (T)/T is extraordi- 
narily large for a transition metal compound (7(1 if) « 
0.42 J/molK 2 ), decreasing rapidly with the temperature 
to w Q.U/molK 2 at 30 K. The spin susceptibility from 
50K to 400K shows a Curie- Weiss [ X = C/{T-6)] behav- 
ior corresponding to weakly antiferromagnetically cou- 
pled (6— -30 to -60 K) vanadium local magnetic moments 
with 5=1/2 and g s» 2, but static magnetic ordering does 
not occur above 0.02 K and superconductivity is not ob- 
served above 0.01 K. The nearly temperature indepen- 
dent spin susceptibility x(T) and Knight shift K(T) for 
T < 30K are a sign of the disappearance of the V local 
moment in this temperature region Q . 

In the traditional heavy fermion compounds there are 
two distinct types of electronic states, the localized f- 
orbitals of lanthanide or actinide atoms which form local 
moments and the delocalized s-, p-, d-orbitals which are 
responsible for the metallic properties. The weak hy- 
bridization of the f-orbitals with the conduction states 
leads to the low temperature anomalies. If one of the 
1.5 3d-electrons per V ion would be in a localized or- 
bital, and the rest in a relatively broad band, then the 
situation would be analogous to f-compounds and all ex- 
perimental facts could be qualitatively understood HQ. 
As the general opinion was that these 1.5 electrons are in 



the (same) band formed by ^g-orbitals, such a model was 
considered to be unrealistic. In this paper we will show 
that the trigonal point group symmetry of the V ion in 
L1V2O4 lifts the degeneracy of the i2g-band*. As a result 
the above mentioned model, is appropriate to estimate 
the heavy-fermion energy scale for this compound. 

The normal-spinel crystal structure is formed by the 
edge-shared oxygen octahedra with V-atoms at the cen- 
tres and Li-atoms between octahedra. The face-centered- 
cubic lattice has four V atoms in the unit cell which form 
a tetrahedron (Fig. |l|). The total space group of the crys- 
tal is cubic but the local point group symmetry of V-ion 
crystallographic position is trigonal. The different trig- 
onal axes of every V-atom in the unit cell are directed 
towards the centre of the tetrahedron. 




FIG. 1. The four V atoms in spinel unit cell with the corre- 
sponding oxygen octahedra. The angular distribution of the 
density of 3d-electrons obtained in the LDA+U calculation 
(mainly the localized A\ g orbitals) is shown. 

The octahedral coordination of the oxygen ions around 
the V results in the strong splitting of d-states into triply 
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degenerate i2g-orbitals with lower energy and double de- 
generate e g -orbitals with higher energy. The band struc- 
ture has three well separated sets of bands, completely 
filled 0-2p-band, partially filled tig band and empty e g 
bands. As only partially filled bands are important for 
the physical properties we will restrict our analysis to the 
t2 g band. 
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FIG. 2. The density of states (DOS) of t 2g band for L1V2O4 
decomposed in Ai g and E g orbitals obtained in the LDA cal- 
culations. 

In the trigonal symmetry crystal field the three ti g 
orbitals are split into the nondegenerate A\ g and double 
degenerate E g representations of the D 3 d group. The 
Ai g -orbital is (xy + xz + yz)/^/3 if the the coordinate 
axes are directed along V-0 bonds. If the z-direction 
is chosen along the trigonal axis then the Aig-orbital is 
3z 2 — r 2 (see Fig. 1) and the E g orbitals have the lobes 
in the plane perpendicular to z. 

We have performed LDA calculations of the electronic 
structure of LiV 2 4 by LMTO method g. In Fig. 2 the 
results for ti g band are presented with the partial DOS 
projected to A\ g and E g orbitals. While the partial E g 
DOS have the width of w 2 eV, the A Xg DOS is much 
narrower (a width of only w 1 eV) . The trigonal splitting 
is smaller than the band width but it is not neglegible. 
We have estimated it as a difference of the centres of 
gravities of A lg and E g DOS's and have found that the 
energy of A\ g orbital is 0.1 eV lower than the energy of 
E g orbitals. 

We have calculated the d-d Coulomb interaction pa- 
rameter U H and have found the value of 3 eV, which is 
larger than the band width W « 2 eV. Such value of the 
U/W ratio leads to the appearance of lower and upper 
Hubbard bands with one electron localized in the for- 
mer and the other 0.5 partially filling the latter. Which 
particular orbital will form the lower Hubbard band is 
defined by the sign of the trigonal splitting between the 



energies of A\ g and E g orbitals. 

The effect of the d-d Coulomb interaction on the elec- 
tronic structure of transition metal compounds can be 
treated by LDA+U method 0. This method is basi- 
cally a combination of the Hartree-Fock approximation to 
multiband Hubbard model with LDA. The ground state 
of the LDA+U solution was found to be indeed a metal 
with one electron nearly completely localized on A\ g or- 
bital and E g orbitals forming relatively broad (« 2 eV) 
band which is partially filled. 

The partial DOS for E g orbitals has a long low-energy 
tail (Fig. 2) and as a consequence the LDA-occupancy 
of Eg orbital is significantly larger than the correspond- 
ing value for Ai g orbital. One would expect that after 
switching on the Coulomb interaction in LDA+U calcu- 
lation, the orbital whose occupation was larger becomes 
in the process of the self-consistency iteration, more and 
more occupied at the expense of all other orbitals and in 
the end the d-electron will be localized on one of the E g 
orbitals. In reality the situation is more complicated. In- 
deed, the Coulomb interaction energy will be lowered by 
localization of the electron on any orbital. However the 
total energy of the solution with the localized electron 
on A\g orbital is lower than energy of the solution with 
Eg orbital due to the trigonal splitting. The rotation 
invariant formulation of the LDA+U method allows the 
system to choose itself on what particular orbital (or the 
particular combination of the basis set orbitals) electrons 
will be localized. If one starts from the LDA orbital oc- 
cupancies, than at the first stage of the self-consistency 
iterations the E g orbital, having larger LDA occupation 
become localized. However further iterations cause 'rota- 
tions' in the five-dimensional space of 3d-orbitals leading 
to the solution with the A\ g orbital occupied. The sys- 
tem arrives at this solution independently of the starting 
point. The total energy as a functional of the orbital 
occupation matrix has only one minimum and the corre- 
sponding LDA+U equations have only one solution. 

The separation of the ^-states into the localized A\ g 
orbital and the conduction band E g orbitals shows that 
LiV204 can be regarded as an analog of f-systems. In 
order to estimate the strength of the interaction between 
the localized and the conduction electrons we have de- 
fined an effective Anderson impurity model. The partial 
DOS for Aig orbital obtained in the LDA+U calculation 
riA lg (E) was used to determine the position of the impu- 
rity state e/ and the hybridization function A(E): 

n Alg (E) = --Zm(E-e f +iA(E))- 1 (1) 

The results are presented on Fig. ||. We then used the 
LDA+U results as input to estimate the Kondo energy 
scale for a single site model. To solve this Anderson 
impurity model we have used a resolvent perturbation 
theory. 
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FIG. 3. The parameters of the effective Anderson impurity 
model. Solid line - hybridization function A(E), dashed arrow 
- the position of the impurity state e/, dashed line - the par- 
tial DOS for Ai g orbital obtained in the LDA+U calculation 
n Alg {E). 

Within this approach one determines the renormal- 
ization of the impurity states by the hybridization with 
the surrounding medium. The latter is given by the hy- 
bridization function A\(E). The approach leads to a per- 
turbation expansion in terms of A(-B). A widely used 
approximation, the Non-Crossing- Approximation (NCA) 
pL is well-known to reproduce the low-energy scale of 
the underlying model which is the single-site Kondo 
temperature Tr-. Here we have extracted this quantity 
directly from the singlet-triplet splitting of the impurity 
states. For the hybridization function A(E) shown in 
Fig. 3 we estimate for the single-site Kondo temperature 
Tk ~ 550K. On the other hand the Kondo temperature 
can be expressed via a Kondo parameter Jk as IlG] : 



(2) 



T K = De 2N (°) j k 



where D is the effective Fermi energy and N(0) the den- 
sity of conduction states per spin. In ref. || the value 
of N(0) was estimated from the experimental value of 
the Pauli susceptibility as 2.9 state/ (eV vanadium spin 
direction). With a D value of approximately 1 eV, this 
gives us an estimation for Jk w 670K. 

In contrast to standard f-systems where the direct ex- 
change between localized f-electrons and conduction elec- 
trons is small, both localized and conduction electrons 
are 3d-electrons in our case and there is a strong onsite 
exchange coupling between them, namely a ferromag- 
netic Hund coupling of the order of leV. The presence 
of two types of d-electrons, namely those in conduction 
states and the electrons forming local moments, is known 
to result in a lattice in a double exchange ferromagnetic 
interaction between local moments. We have estimated 
this value of the intersite exchange coupling parameter 
from the results of LDA+U calculation using the formula 
derived as a second derivative of the total energy in re- 
spect to the angle between local moments directions [[llj . 
Our calculation gave the value of the double exchange 
parameter Jdex= 530K. 



However, there is another contribution to the interac- 
tion between the local moments. The Kondo exchange 
is between a local moment (electron on the Ai 9 -orbital) 
on one site and the spin of the conduction electron (E g 
-orbital) on the neighboring site, because having differ- 
ent symmetry the orbitals do not mix on the same site. 
As the spins of the E g and A\ g electrons on the neigh- 
boring site are 'strongly' coupled by the Hund interac- 
tion, it gives us an effective antiferromagnetic (AF) in- 
teraction between local moments on neighboring sites ap- 
proximately equal to the Kondo exchange parameter Jk- 
As a result there is a strong cancellation between these 
two processes, ferromagnetic double exchange and AF 
Kondo-induced exchange so that the net exchange inter- 
action is small. Simply substracting the two terms leads 
to an estimate Jx-Jdex ~ 140K. The measured Curie- 
Weiss O at high temperature (O ~ 5QK) points to a 
more complete cancellation and a net AF exchange in- 
teraction that is an order in magnitude smaller than our 
estimate. This discrepancy is not surprising in view of 
the difficulty of making a reliable estimate in the presence 
of this strong cancellation. Note the inherent frustration 
that inhibits the onset of AF-order in a spinel lattice 
does not enter the determination of the value of O. The 
small value of O shows that the net exchange interaction 
between neighboring local moments is very weak. 

A realistic lattice model for LiV204 contains two com- 
peting terms which couple the conduction and localized 
states. These are the onsite Hund's ferromagnetic cou- 
pling and the AF Kondo interaction which couples con- 
duction and localized electrons on neighboring sites. This 
competition makes difficulties for a first principles treat- 
ment. On the other hand, there are several arguments 
in favor of ignoring the ferromagnetic interactions rela- 
tive to the Kondo interaction. For a single localized site, 
it is well-known that the onsite ferromagnetic coupling 
between conduction and localized states scales to weak 
coupling limit at low temperatures and so can be ignored. 
However it could be argued that in the lattice this fer- 
romagnetic coupling scales to the strong coupling limit 
through the well-known double exchange effect. But as 
discussed above the double exchange effect is cancelled 
here by the AF interaction between the localized spins 
induced by the Kondo effect. Therefore it seems plau- 
sible to ignore at least as a first step the ferromagnetic 
interactions and treat only the AF interaction so that the 
model is simply a Kondo lattice model. There may be 
some renormalization of the Kondo exchange parameter, 
Jk but for now we ignore that too. 

In the Kondo lattice model there is also a competi- 
tion between the induced AF interactions which favor 
AF order and the Kondo effect which favors a singlet 
groundstate |l2|. Here, the former are very weak due to 
the cancellation affects discussed above and as a result 
we are in the limit where the Kondo effect dominates. 
This leads to the formation of a heavy-fermion Landau 
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Fermi liquid with a characteristic temperature scale for 
the onset of quantum coherence, T co h- The exact value 
of T co h is difficult to estimate but there are a number of 
strong arguments that T co h <C Tk, the single site Kondo 
temperature. We summarize these arguments below. 

In the single site case, the local moment forms a sin- 
glet pair with any conduction electron within an energy 
Tk of the Fermi energy. We have the picture of a com- 
plex screening cloud which delocalizes the impurity at 
the Fermi level. The resulting non-perturbative ground 
state is found to be of Fermi-liquid type |l3|,[l4|], where 
all the physical quantities depend on Tk- In the con- 
centrated Kondo lattice, the number of conduction elec- 
trons per site to screen the local moments is of order 
nk = -jj- -C 1. There is a lack of conduction electrons 
to screen all the spin array at the energy scale Tk- this 
is the well-known exhaustion phenomenon p5| . In that 
sense, a macroscopic singlet ground state should not take 
place at Tk, but rather at a very low temperature T co h- 
To obtain it, we can use the following simple thermo- 
dynamical arguments. The condensation energy that we 
really dispose to screen the (effective) localized spins is: 
E con d = NynfeT/f where Ny is the total number of V- 
ions (we have one localized 3d-electron per V) . In view to 
stabilize a perfect singlet ground state, this energy must 
absorb all the entropy of the impurity lattice which is de- 
fined as: Stot — Ny In 2. Then, the energy scale at which 
the entropy of the spin array will go to zero can be simply 
defined as T coh = E cond /S tot ~ n k T K = T K 2 / D |§. For 
T <C T co h, the prevalent bonds should be formed between 
local moments and then T co h plays the role of the effec- 
tive Kondo temperature in the lattice problem jl7| . It 
should be noted that these effects predicted by physical 
arguments could recently be indeed observed in calcu- 
lations for the periodic Anderson model in the frame- 
work of the Dynamical Mean-Field Theory [H . Taking 
into account the estimations for Tk and D, we obtain 
T co h ~ 25 — 4QK, which agrees well with the numerical 
result for the periodic Anderson model |n| . It also seems 
to be in good agreement with the experimental result in 
L1V2O4. This should produce an enhanced linear spe- 
cific heat and Pauli susceptibility which are proportional 
to Ny /T C oh, and a (normalized) Wilson ratio Rw which 
is equal to one |l7| . Experimentally, Rw is on the order 
of unity as in conventional heavy- metals p| . 

In conclusion, our calculations using the LDA+U 
method give a theoretical justification of a model with 
one of the 3d-electrons per V localized in a ^4i s -orbital 
and the remaining 0.5 electron/V in a conduction band 
state, primarily of Eg symmetry, which has been pre- 
viously discussed ||l4|]. This leads to a lattice model 
with competing onsitc ferromagnetic coupling due to the 
Hund's rule and nearest neighbor antiferromagnetic cou- 
pling due to the Kondo effect. We present arguments 
that such a model reduces to a Kondo lattice model. Es- 
timates for the temperature scale for the onset of quan- 



tum coherence give a small value, much less than the 
single site Kondo temperature, in agreement with exper- 
iments. The low temperature heavy fermion Fermi liquid 
is strongly correlated and is therefore a good candidate 
for a transition to unconventional superconductivity, if it 
can be made perfect enough. 

*Note added — Recently an LDA calculation of the elec- 
tronic structure of LiV2C>4 was made by Eyert et al. |2C| ] 
where the partial density of states of t2 g band was ana- 
lyzed using trigonal symmetry but the possibility of or- 
bital polarization of electrons due to the Coulomb inter- 
action was not investigated. 
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